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Abstract 

If quantum chromodynamics (QCD) is renormalized by minimal subtraction 
(MS), at higher orders, the strong couphng constant Og and the quark masses ruq 
exhibit discontinuities at the flavour thresholds, which are controlled by so-called 
decoupling constants, and (m, respectively. Adopting the modifled MS (MS) 
scheme, we derive simple formulae which reduce the calculation of and to the 
solution of vacuum integrals. This allows us to evaluate (g and Cm through three 
loops. We also establish low-energy theorems, valid to all orders, which relate the 
effective couplings of the Higgs boson to gluons and light quarks, due to the virtual 
presence of a heavy quark h, to the logarithmic derivatives w.r.t. nih of Cg and Cm, 
respectively. Fully exploiting present knowledge of the anomalous dimensions of 
and TUq, we thus calculate these effective couplings through four loops. Finally, we 
perform a similar analysis for the coupling of the Higgs boson to photons. 
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1 Introduction 



As is well known, in renormalization schemes based on the method of minimal subtraction 
(MS) IjlJ, including the modified MS (MS) scheme 0, which is routinely used in quantum 
chromodynamics (QCD), the Appelquist-Carazzone decoupling theorem does not hold 
true in its naive sense. Let us consider QCD with ni = nj — 1 light quarks flavours q and 
one heavy flavour h. Then, the contributions of h to the Green functions of the gluons 
and light quarks expressed in terms of the renormalized parameters of the full theory do 
not exhibit the expected 1/m/j suppression, where rrih is the mass of h. The reason for 
this is that the (3 and 7^ functions governing the running of the strong coupling constant 
as{fi) and the light-quark masses mq{fi) with the renormalization scale fi do not depend 
on any mass. 

The standard procedure to circumvent this problem is to render decoupling explicit by 
using the language of effective fleld theory, i.e. h is integrated out. Speciflcally, one con- 
structs an effective ra^-flavour theory by requiring consistency with the full rij-flavour the- 
ory at the heavy-quark threshold /i*^"^^ = 0{mh) @J^. This leads to nontrivial matching 
conditions between the couplings and light-quark masses of the two theories. Although, 
Q,(n;)(i7^^) = ai^-^\mh) and m^^^\mh) = rn~q^\mh) at leading (tree-level) and next-to- 
leading (one-loop) order, these identities do not generally hold at higher orders in the MS 
scheme. Starting at next-to-next-to-leading (two-loop) order, they are broken by flnite 
corrections, of C(af), as was noticed in the pioneering works of Refs. @,|]. The relations 
between the couplings and light-quark masses of the full and effective theories are called 
decoupling relations; the proportionality constants that appear in these relations are de- 
noted decoupling constants, C3 and Cm- In this paper, and Cm are computed through 
next-to-next-to-next-to-leading (three-loop) order, 0{a^g). They have to be applied when- 
ever a flavour threshold is to be crossed. If the /i evolutions of ai"^^(/x) and mg"-^^(/i) are 
to be performed at loops, then consistency requires that the decoupling relations be 
implemented at A^ — 1 loops. Then, the residual /i dependence of physical observables will 
be of the (A^ -|- l)-loop order. If our new results are combined with the recently evaluated 
four-loop coefficients of the (3 and 7^ 101 functions, then it is possible to consistently 
describe QCD-related observables through 0{a^) 

The dominant subprocess for the production of the standard-model (SM) Higgs boson 
at the CERN Large Hadron Collider (LHC) will be the one via gluon fusion. Therefore, 
an important ingredient for the Higgs-boson search will be the effective coupling of the 
Higgs boson to gluons, usually called Ci. A standard approach to study this coupling as 
well as the Higgs-boson coupling to photons has been to use low-energy theorems (LET's) 
0. At the two- loop level, Ci has been known since long |jlO|,|Tl[. Recently, the three- loop 
correction to Ci has been obtained through a direct diagrammatic calculation [O. In 



Ref. |]I2[, also a LET which allows for the computation of Ci from the knowledge of C,g has 
been used. In a similar way, the effective couplings of the Higgs boson to light quarks may 
be treated as well. The corresponding coefficient function, which is usually denoted by 
C21 has a similar connection to Cm- In this paper, we establish the relationships between 
Ci,C2 and CgXm to all orders by formulating appropriate LET's. With their help, we 
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compute Ci and C2 through four loops, i.e. 

The outhne of the paper is as follows. In Section |], we derive simple formulae which 
allow us to determine the decoupling constants and Cm for ^^(/i) and respectively, 
by simply evaluating vacuum integrals. With the help of these formulae, in Section |^, we 
calculate and Cm up to the three-loop order. In Section ^ we derive LET's, valid to all 
orders, for the coefficient functions Ci and C2, which comprise the virtual top-quark effects 
on the interactions of the Higgs boson with gluons and quarks, respectively. Exploiting 
these let's, we compute Ci and C2 through four loops. In Section ^ we extend our 
analysis to also include the interaction of the Higgs boson with photons. In Section |^, we 
explore the phenomenological significance of our results and present our conclusions. In 
Appendix 0, we display the dependence of our key results, adopting the general gauge 
group SU(Ai'c). In Appendix ^, we list the decoupling constants for the quark and gluon 
fields in the covariant gauge. 



2 Formalism 



Let us, in a first step, fix the notation and present the framework for our calculation. 
Throughout this paper, we work in the MS renormahzation scheme 0. In the main part 
of the paper, we concentrate on the QCD gauge group SU(3), i.e. we put Nc = 3. The /3 
and 7m functions of QCD are defined through 



(«/) 



fi^d as 
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where denotes the number of loops. Recently, the four-loop coefficients /^a"'-'^^ and 7^,3'' 
have become available MM- According to present knowledge, we have 
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where ( is Riemann's zeta function, with values ({2) = 7r^/6, C(3) ~ 1.202 057, ({4) 
7rV90, and C(5) ^ 1.036 928. 

The relations between the bare and renormalized quantities read 
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where Qs = \/47m^ is the QCD gauge coupling, /x is the renormahzation scale, D = A — 2e 
is the dimensionality of space time, ipq is a quark field with mass m^, GJJ is the gluon 
field, and c" is the Faddeev-Popov-ghost field. For simplicity, we do not display the 
colour indices of the quark fields. The gauge parameter, ^, is defined through the gluon 
propagator in lowest order. 



q-^ + ie 
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The index '0' marks bare quantities. We should mention that relations @ hold true both 
in the full n/-flavour and effective nrflavour theories. Only the renormahzation constants 
Zg and Zm are relevant for our purposes. They are known through order 0,01 • The 

other relations are only listed in order to fix the notation. 

In a similar way, we can write down the relations between the parameters and fields 
of the full and effective theories. From the technical point of view, it is simpler to first 
consider these relations for bare quantities and to construct those for the renormalized ones 
afterwards. The bare decoupling constants are defined through the following equations: 
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where the primes mark the quantities of the effective n^-flavour theory. 

Combining Eqs. and (^, we obtain relations between the renormalized coupling 
constants and masses of the full and effective theories, viz 



(7) 
(8) 
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Here, it is understood that the right-hand sides are expressed in terms of the parameters 
of the full theory. The unknown quantities in these decoupling relations are (g and Cm- 

In a next step, let us introduce the effective Lagrangian, C, which depends on the 
decoupling constants and on the bare parameters and light fields of the full theory. It 
is obvious from gauge invariance that the most general form of C is simply the one of 
usual QCD, C^'^^, retaining, however, only the light degrees of freedom. Specifically, the 
definition of C reads 

^' {g's, ml e-, c°'- c°) = {g'j. mi', e-, <, cj'^ c°'"^) , (9) 

where q represents the rii light quark flavours and collectively denotes all bare decou- 
pling constants of Eq. (^. Exploiting the circumstance that the result for some n-particle 
Green function of light fields obtained from £Q'^° in the full theory agrees, up to terms 
suppressed by inverse powers of the heavy-quark mass, with the corresponding evaluation 
from C in the effective theory, we may derive relations which allow us to determine the 
decoupling constants (f. 

As an example, let us consider the massless-quark propagator. Up to terms of (9(1/ rrih), 
we have 

1 



(10) 
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^[i + so;(p2)]' 

where the subscript 'V reminds us that the QCD self-energy of a massless quark only 
consists of a vector part. Note that Sy(p^) only contains light degrees of freedom, whereas 
Ey(p^) also receives virtual contributions from the heavy quark h. As we are interested in 
the limit rrih — oo, we may nullify the external momentum p, which entails an enormous 
technical simplification because then only tadpole integrals have to be considered |l^ . 
Within dimensional regularization, one also has Sy(0) = 0. Thus, we obtain 

C° = i + s°^(o), (11) 

where the superscript h indicates that only the hard part of the respective quantities needs 
to be computed, i.e. only the diagrams involving the heavy quark contribute. 

In a similar fashion, it is possible to derive formulae for the decoupling constants 
(^3, and as well. These read 
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C° = i + nf(o), (14) 



5 



where Sy(p^) and ^^(p^) are the vector and scalar components of the hght-quark self- 
energy, defined through E(p) = ^Ey(p^) +mgE5(p^), and Ucip^) and Hdp^) are the gluon 
and ghost vacuum polarizations, respectively. Specifically, Ucip"^) and Hdp^) are related 
to the gluon and ghost propagators through 
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ab 
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dxe'P-'' (Tc''^\x)(P^\Q)), (15) 



respectively. Finally, we extract an expression for C^^ from the relationship between the 
bare Gcc-vertex form factors of the full and effective theories. 
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by nullifying the external four-momenta p and k. Here, r^g^(p, k) is defined through the 
one-particle-irreducible (IPI) part of the amputated Gcc Green function as 
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where p and k are the outgoing four-momenta of c and G, respectively, and /"''^ are the 
structure constants of the QCD gauge group. We thus obtain 
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where 



C? = i + r°f,,(o,o). (19) 

In contrast to the renormalization constants Zi in Eq. (1), the decoupling constants 
also receive contributions from the finite parts of the loop integrals. Thus, at C(a^), 
we are led to evaluate three-loop tadpole integrals also retaining their finite parts. 



3 Decoupling relations 

In this section, we compute the bare decoupling constants and Cm combine them 
with the well-known 0{a1) results for Zg and Zm to obtain the finite quantities Cg and 
Cm- As may be seen from Eq. (p!2D , the computation of Cm requires the knowledge of the 
hard part of the light-quark propagator. There are no one-loop diagrams contributing to 
Sy'(O) and TPg{Q). At two loops, there is just one diagram. Even at three loops, there 
is only a moderate number of diagrams, namely 25. Typical specimen are depicted in 
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EfiO). It 



Fig. |I]. Actually, through three loops, Eq. (JT^) simplifies to Cm = 1 
should be noted that the vector and scalar parts separately still depend on the QCD gauge 
parameter ^, but ^ drops out in their sum, which is a useful check for our calculation. We 
use the program QGRAF [1^ to generate the relevant diagrams. The computation is then 



performed with the help of the program MATAD ||T5|], which is based on the technology 
developed in Ref. [|16| and written in FORM Using Eq. (P), we finally obtain 
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with Li4 being the quadrilogarithm, is a constant typical for three-loop tadpole diagrams, 
rii = Uf — 1 is the number of light-quark fiavours q, and rrih = Tn^f^^\n) is the MS mass 
of the heavy quark h. For the numerical evaluation in the last line of Eq. (PU]), we have 

chosen /i = /i/j, where fih = '^^h^\l^h)- The C(af) term of Eq. (|20|) agrees with Ref. 
the 0{a1) term represents a new result. The generalization of Eq. ( PP]) appropriate for 
the gauge group SU(Ai'c) is listed in Appendix 0. 

For the convenience of those readers who prefer to fix the matching scale /i in units 



of the pole mass, M^, we substitute in Eq. (^) the well-known relation between 
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Figure 1: Typical three-loop diagrams pertinent to Sy'(O) and T,g^{0). Solid, bold-faced, 
and loopy lines represent massless quarks q, heavy quarks h, and gluons G, respectively. 



According to Eq. ([TsD, three ingredients enter the calculation of Cg, namely the hard 
parts of the gluon and ghost propagators and the gluon-ghost vertex correction. At one 
loop, only one diagram contributes, namely the diagram where the gluon splits into a 
virtual pair of heavy quarks. At two loops, three diagrams contribute to 11^'* (0) and 
one to n|?^(0). In the case of r^^(0,0), there are five diagrams, which, however, add 
up to zero. To this order, the three contributions are still separately independent of the 
gauge parameter ^, so that the ^ independence of their combination does not provide a 
meaningful check for our calculation. The situation changes at 0{al), where all three parts 
separately depend on ^ and only their proper combination is ^ independent as is required 
for a physical quantity. At this order, the numbers of diagrams pertinent to n|^(0), 
ll^^{0), and r^g^(0,0) are 189, 25, and 228, respectively. Typical representatives are 
shown in Fig. The complexity of the problem at hand necessitates the use of powerful 
analytic technology |]T^^|I6],|1^] to organize the calculation. Inserting into Eq. (|^) the 
result for thus obtained finally leads to the following answer in the MS scheme: 
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The 0{al) term in Eq. (^3|) has recently been published |P. Leaving aside this term, 
Eq. (p3D agrees with Ref. [|19[, while the constant term in 0{a1) slightly differs from the 
result published in Ref. In the meantime, the authors of Ref. P] have revised 
their original analysis and have found agreement with Ref. |19[. The SU(A^c) version of 
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faced, loopy, and dashed lines represent heavy quarks h, gluons G, and Faddeev-Popov 
ghosts c, respectively. 



Eq. ( p3D may be found in Appendix |^. Introducing the pole mass Mh leads to 
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Notice that the 0{al) terms of Cm and Q depend on the number n; of light (massless) 
quark flavours. However, this dependence is feeble. 

In the framework of the QCD-improved parton model, quarks and gluons appear as 
external particles, so that the knowledge of the decoupling constants C2 and Cs, which 
emerge as by-products of our analysis, is actually of practical interest for higher-order 
calculations. In contrast to Cm and Cg, C2 and Cs are ^ dependent. For future applications, 
we list C,2 and Cs in Appendix |B[ 



4 Low-energy theorems for the ggH and qqH inter- 
actions 

The interactions of the SM Higgs boson with gluons and light quarks are greatly affected 
by the virtual presence of the top quark. In fact, the Higgs-boson coupling to gluons is 
essentially generated by a top-quark loop alone. In general, the theoretical description of 
such interactions is very complicated because there are two different mass scales involved. 
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Mh and Mt. However, in the limit Mh <^ 2Mt, the situation may be greatly simplified by 
integrating out the top quark, i.e. by constructing a heavy-top-quark effective Lagrangian. 

The starting point of our consideration is the bare Yukawa Lagrangian of the full 
theory, 

^0 "/ 



(25) 



which governs the interactions of the neutral CP-even Higgs boson H with all nj quark 
fiavours, including the heavy one. Here, v is the Higgs vacuum expectation value. The 
heavy-quark effective Lagrangian describing the interactions of H with the gluon G and 
the rii light-quark fiavours may be written in bare form as 
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The operators 0[ are only constructed from light degrees of freedom and read [0j2l|j22 
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where G^^, is the colour field strength, = — 'i/x^^(7s(A'^/2)G^ and V^* = 5°-^d^ 



ilJ?^gsf°'^'^G'j^ are the covariant derivatives acting on the quark and gluon/ghost fields, 
respectively, and A° are the Gell-Mann matrices. The residual dependence on the mass 
iTLh of the heavy quark h is contained in the coefficient functions C°. 

In phenomenological applications, one is mainly interested in the renormalized co- 
efficient functions Gi and C2, since only these contribute to physical observables. The 
renormalization of and 0[ (z = 1,2) has been explained in Ref. [^. For the reader's 
convenience, we repeat here the key results. Denoting the renormalized operators by 
square brackets, we have ETI, 
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Note that Og, O'^, and O'^ do not mix with O'^ and 0'2- On the other hand, the coefficient 
functions are renormalized according to 
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Consequently, the physical part of £ofr takes the form 



C^^' = -^{C,[0[]^C,[0',]). (31) 



Ci and [0'j\ {i = 1, 2) are individually finite, but, with the exception of [02], they are not 



separately renormalizat ion-group (RG) invariant. In Ref. [^, a RG-improved version of 
Eq. ([31| ) has been constructed by exploiting the RG-invariance of the trace of the energy- 
momentum tensor. The ratio H^/v^ receives a finite renormalization factor, which is of 
0{GfM^). Its two- and three-loop QCD corrections have been found in Refs. |^ and 
respectively. 



The derivation of formulae to compute the coefficient functions is very similar to 
the case of the decoupling relations. As an example, let us consider the derivation of a 
formula involving C° and C^. The starting point is the IPI Green function of two gluons 
which contains a zero-momentum insertion of the composite operator Oh = m^lhPh^ . In 
momentum space, it reads in bare form 



GGO; 
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g'^'^P^^GGoSP'^) + terms proportional to p^p'^ , (32) 



where p denotes the four-momentum flowing along the gluon line. In the limit rrih 00, 
Oh niay be written as a linear combination of the effective operators given in Eq. (P?]), so 
that 

^GGO- ip) = t f dxdye'P<--y^ (TG'''''>^ix)G'''-''iy) (C^O\ + CM)^' + 



= ^Cs / dxdy e'P-^^'-y^ (TG°''"'^(x)G°'''^'^(y) [G°,0[ + G^O'^))^^^ + ... 

= -^''VC3(-4C° + 2C°) (1 + higher orders) + . . . , (33) 

where the ellipses indicate terms of 0{l/mh) and terms proportional to p'^p". In the 
second step, we have used Eq. (^ together with the fact that ^c'cOhiP) represents an 
amputated Green function. If we consider the coefficients of the transversal part in the 
limit p — 0, we observe that the contributions due to the higher-order corrections on the 
r.h.s. of Eq. (|33| ) vanish, as massless tadpoles are set to zero in dimensional regularization. 
The contribution due to the other operators vanish for the same reason. On the l.h.s., 
only those diagrams survive which contain at least one heavy-quark line. Consequently, 
the hard part of the amputated Green function is given by 

r^GO,(0) = C3°(-4C^? + 2C°). (34) 
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It is convenient to generate the diagrams contributing to the l.h.s. of Eq. (0) by dif- 
ferentiating the gluon propagator w.r.t. the heavy-quark mass, so that we finally arrive 
at 

C3°(-4C? + 2C°) = -Id'.U'iiO), (35) 
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where 9° = [(m°)^9/c}(m° 

coupling of the Higgs boson to the heavy quark is proportional to m° 



The last equation results from the fact that the Yukawa 
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In a similar way, we obtain four more relationships, namely 
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We may now solve Eqs. (pSj) and (^) for the coefficient functions C° (i = 1, . . . , 5). 

It is tempting to insert Eqs. (|TT])-(0) into Eqs. ( P5| ) and (|5BD, so as to express the 
bare coefficient functions in terms of derivatives of the decoupling constants w.r.t. m°. 
Solving the three equations involving C°, C^, and for C°, we obtain 
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(37) 



Next, we express C° through renormalized quantities. Using d'^ = dh, we find 
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2 



c^hlna'. 



(38) 



Identifying the renormalization factor of Eq. 



we obtain the amazingly simple relation 



(39) 



This relation opens the possibility to compute C\ through (9(a^), since one only needs 
to know the logarithmic pieces of in this order, which may be reconstructed from its 
lower-order terms in combination with the four-loop /3 and 7^ [0 functions. It is 
possible to directly relate C\ to the /3 and 7^ functions of the full and effective theories. 
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Exploiting the relation 
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where a'^ = a'g{fi, a^, m/j), we find 
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In the case of C2, we may proceed along the same lines to obtain 
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where m'^ = m'g{fi,as,mh). It should be stressed that Eqs. (^If ) and (^2]) are valid to all 
orders in a^. 

Fully exploiting present knowledge of the /5 and 7^ [0 functions, we may evaluate 
Ci and C2 through 0{a'l) via Eqs. (^T]) and (|42|). In the pure MS scheme, we find 
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where, for simplicity, we have chosen /x = /i/j in the approximate expressions. The 0{al) 
term of Eq. (i4) may also be found in Ref. |26]. The corresponding expressions written 
with the pole mass Mh read 
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where we have put fi = Mh in the numerical evaluations. The 0{al) term of Eq. 
may also be found in Ref. |]12|, where the hadronic decay width of the SM Higgs boson 
has been calculated through C(a^). 



5 Low-energy theorem for the 77// interaction 

In this section, we extend the formalism developed in Sections ^ and ^ to include the 
Higgs-boson interactions with photons. From the technical point of view, this means 
that we are now concerned with QCD corrections to Green functions where the external 
particles are photons instead of gluons. In the following, this is indicated by an additional 
subscript '7'. Furthermore, we do not need to consider Green functions involving external 
ghost lines any more. The Abelian versions of Eqs. (^), (^, and (^) emerge via the 
substitutions g ^ e = y/InE, G"^ A^, F^^, jl 1, and f"-'"' 0, where e is 

the gauge coupling of quantum electrodynamics (QED), d is the fine-structure constant, 
Afj^ is the photon field, and F^^, is the electromagnetic field strength. We continue to work 
in the MS renormalization scheme. The system of composite operators ([27| ) is reduced 
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to the Abelian counterparts of the gauge- invariant operators O'l, O2, and O'^. The other 
formulae derived in Sections |^ and § simphfy accordingly. 

Similarly to Eq. (|^), the decoupling relation for the renormalized MS fine-structure 
constant reads 
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To three loops in QCD, the photon wave-function-renormalization constant Z^^ may be 
extracted from Ref. |2^, while the bare decoupling constant (^^ for Q;(/i) is determined by 
the hard part of the three-loop photon propagator, which may be found in Refs. p8|j29| . 
Putting everything together, we find in the pure MS scheme 
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where Qg is the fractional electric charge of quark flavour q. Introducing the pole mass 
M/i, this becomes 
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Next, we turn to the coefficient function Ci^ of the heavy-quark effective 77if couphng. 
Equation (^l]) undergoes obvious modifications to become 
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where f3^ is the (3 function governing the running of the MS fine-structure constant a{fi). 
An expression for (3^ may be extracted from the well-known QCD corrections to the 
photon propagator 0]. Restricting ourselves to the leading order in a, we find 
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We are now in the position to evaluate Ci^ up to the four-loop level. In the pure MS 
scheme, we find 
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The corresponding expression written with the pole mass Mh reads 
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The C(q:Q!s) term of Eq. (|56|) agrees with Ref. ||3T|. By comparing Eqs. ( ^5]) and (p6D, we 
notice that only the 0{aal) terms depend on the definition of the top quark mass. The 
0{aa1) contributions due to the diagrams where the Higgs boson and the photons are 
connected to the same quark loop agree with Ref. where the corresponding vertex 
diagrams have been directly calculated. In Ref. also power-suppressed terms of 

0{aa1) have been considered. 



6 Discussion and conclusions 

In the context of QCD with MS renormalization, the decoupling constants Cm, and C,g given 
in Eqs. (^) and (0), respectively, determine the shifts in the light-quark masses mq{fi) 
and the strong coupling constant as{fi) that occur as the threshold of a heavy-quark 
fiavour h is crossed. In Eqs. (|20| ) and (p3D, the matching scale /i, at which the crossing 
is implemented, is measured in units of the MS mass m/j(/i) of h. The corresponding 
decoupling relations, formulated for the pole mass M/j, are given in Eqs. ( [2^ ) and (P^. 
Equations (|20|) and (|^) are valid through three loops, i.e. they extend the results of 
Refs. [^|l9l by one order. These results will be indispensible in order to relate the QCD 
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predictions for different observables at next-to-next-to-next-to-leading order. Meaningful 



estimates of such corrections already exist [p^] . 

We wish to stress that our calculation is based on a conceptually new approach that 
directly links the decoupling constants to massive tadpole integrals. This is crucial because 
the presently available analytic technology does not permit the extension of the methods 
employed in Refs. to the order under consideration here. In fact, this would require 

the evaluation of three-loop diagrams with nonvanishing external momenta in the case of 
Ref . , or even four-loop diagrams in the case of Ref . ]19[ . 

We are now in a position to explore the phenomenological implications of our results. 
To that end, it is convenient to have perturbative solutions of Eqs. (|l|) and (0) for fixed 
Uf in closed form. Iteratively solving Eq. (0) yields IQ 
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where bjy = Pn/Po {N = 1,2,3), L = ln(/i^/A^), and terms of C(l/L^) have been ne- 
glected. The asymptotic scale parameter. A, is defined in the canonical way, by demanding 
that Eq. ( |57| ) does not contain a term proportional to (const. /L^) Equation ( p7D ex- 
tends Eq. (9.5a) of Ref. |3^ to four loops. Combining Eqs. (P and (H), we obtain a 
differential equation for mg{^) as a function of as(/i). It has the solution UTI 



mg(/io) c(as(/io)/vr) 
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(59) 



where cat = Jn/Po = 0, . . . , 3) and terms of O^x"^) have been neglected. 

Going to higher orders, one expects, on general grounds, that the relation between 



(■n-l)( 



{jj,') and ah'^\iJt)^ where /x' -C /i^"-''^ <^ jj. 



becomes insensitive to the choice of the 
. This has been checked in Ref. for 



matching scale, as long as /i*^"'-^^ = 0{mh 

three-loop evolution in connection with two-loop matching. Armed with our new results, 
we are in a position to explore the situation at the next order. As an example, we consider 
the crossing of the bottom-quark threshold. In particular, we wish to study how the jj^^^ 
dependence of the relation between a^'^^Mr) and a^f\Mz) is reduced as we implement 
four-loop evolution with three-loop matching. Our procedure is as follows. We first 
calculate a^J^\jj^^^) by exactly integrating Eq. (|l|) with the initial condition a^^^M^-) = 
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af(M,) 




(5) 



Figure 3: /i^^^ dependence oia'f^Mz) calculated from a'f^Mr) = 0.36 and Mb = 4.7 GeV 
using Eq. (|l|) at one (dotted), two (dashed), three (dot-dashed), and four (solid) loops in 
connection with Eqs. (|^ and (|2^ ) at the respective orders. 



0.36 0, then obtain af)(/x(5)) from Eqs. (0) and (^ with Mb = 4.7 GeV, and finally 
compute af'^Mz) with Eq. (|l|). For consistency, A^-loop evolution must be accompanied 
by (A^ — l)-loop matching, i.e. if we omit terms of 0{af'^'^) on the right-hand side of 
Eq. ([1|), we need to discard those of 0{a^) in Eq. (^) at the same time. In Fig. ^, the 
variation of a^^^Mz) with fi^^^ /Mb is displayed for the various levels of accuracy, ranging 
from one-loop to four-loop evolution. For illustration, fi^^'' is varied rather extremely, 
by almost two orders of magnitude. While the leading-order result exhibits a strong 
logarithmic behaviour, the analysis is gradually getting more stable as we go to higher 
orders. The four-loop curve is almost flat for /i^^^ ^ 1 GeV. Besides the n^^^ dependence 
of a^^\Mz), also its absolute normalization is significantly affected by the higher orders. 
At the central matching scale /i^^-* = Mb, we encounter a rapid, monotonic convergence 
behaviour. 
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mf^ (M J 




(5) 



Figure 4: fi^^^ dependence of m!^^\Mz) calculated from /x, 



(5), 



m 



(4) 



(/ie) = 1.2 GeV, Mb 



4.7 GeV, and af\Mz) = 0.118 using Eq. (|9D at one (dotted), two (dashed), three (dot- 
dashed), and four (solid) loops in connection with Eqs. (||) and ( p2D at the respective 
orders. 



Similar analyses may be performed for the light-quark masses as well. For illustration, 
let us investigate how the n^^^ dependence of the relation between /ic = m^^^Hc) and 
mf^Mz) changes under the inclusion of higher orders in evolution and matching. As 
typical input parameters, we choose /ic = 1.2 GeV, Mi, = 4.7 GeV, and a^^\Mz) = 0.118. 
We first evolve m^^\fi) from = /ic to /i = fi^^^ via Eq. (^), then obtain m^^\fi^^^) from 
Eqs. (H) and 
„("/), 



2]), and finally evolve -'(/i) from /i = fi^ ' to /i 



Mz via Eq. (H). In all 

steps, as'^'ifJ') is evaluated with the same values of nj and /x as nic'^\ii). In Fig. H, we 
show the resulting values of m^^)(M^) corresponding to A^-loop evolution with (A^ — 1)- 
loop matching for N = 1,...,4. Similarly to Fig. ^ we observe a rapid, monotonic 
convergence behaviour at the central matching scale /i'-^^ = Mf,. Again, the prediction for 
A^ = 4 is remarkably stable under the variation of /i'-^-* as long as /i'-^-* ^ 1 GeV. 

An interesting and perhaps even surprising aspect of the decoupling constants (g and 
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Cm is that they carry the full information about the virtual heavy-quark effects on the 
Higgs-boson couplings to gluons and light quarks, respectively. In fact, the relevant 
coefficient functions Ci and C2 in the renormalized version of the effective Lagrangian (p6|) 
emerge from (g and Cm, respectively, through logarithmic differentiation w.r.t. the heavy- 
quark mass ruh- In the LET's (^) and (0), these relationships have been established 
to all orders. By virtue of these LET's, we have succeeded in obtaining the four-loop 
corrections to Ci and C2 from the three-loop 0{a^) expressions for (g and Cm 
complemented by the four-loop results for the P and 7m functions. The pure 
MS expressions for Ci and C2 may be found in Eqs. (^31) and (0), respectively; the 
corresponding versions written in terms of the heavy-quark pole mass are listed in 
Eqs. dl) and (H). 

From the phenomenological point of view, Ci is particularly important, since it enters 
the theoretical prediction for the cross section of the gg ^ H subprocess, which is expected 
to be the dominant production mechanism of the SM Higgs boson H at the CERN proton- 
proton collider LHC. The C(a^) term of Eq. ( ^3]) is routinely included in next-to-leading- 
order calculations within the parton model Recently, a first step towards the next- 
to-next-to-leading order has been taken in Refs. |3^j3^ by considering the resummation 
of soft-gluon radiation in pp H + X. The 0{al) formula (13) for Kh^n in the revised 
version of Ref. agrees with our Eq. (^31) for c[^^\jj,) to this order, if we identify 
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1 «f ) 



(Mh) 
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TT 



(60) 



The RG improvement of C\ in the first line of Eq. (^) is adopted from Eq. (30) of 
Ref. p3[, where it has been employed in the context of three-loop 0{p?fipMl^ corrections 
to hadronic Higgs-boson decays. We stress that the particular choice of scales in Eq. (|60|) 



is essential in order to recover Eq. (13) of Ref. |36| from the general framework elaborated 
here. The /i dependence of k^^h is not displayed in Ref. We note that the O(a^) 

result for C\ was originally obtained in Ref. |1^ . 



In their Eq. (8), the authors of Ref. |^ also present an alternative version of a generic 
formula for C\\ see also Eq. (25) of Ref. A central ingredient of this formula is 

the quantity /3t(as), which they paraphrase as the top-quark contribution to the QCD (5 
function at vanishing momentum transfer. According to Ref. [^, i^tip-s) may be obtained 
from some gauge-independent variant of the top-quark contribution to the bare gluon self- 
energy at vanishing momentum transfer. Since we are unable to locate in the literature, 
including Refs. |]36| , |37|1 and the papers cited therein, a proper general definition of /9t(as) 
in terms of the familiar quantities of QCD, it remains unclear whether Eq. (8) of Ref. [0^ 
agrees with our Eq. (^) in higher orders. If we were to fix fitip-s) so that Eq. (8) of 
Ref. I3H] reproduces our Eq. (ffTl) after the RG-improvement of Eq. (|60D, the result would 
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be 



where a[^^ 



(61) 



m 



(6) 



We emphasize that, in order to derive Eq. (PT|), we found 
it indispensable to consider the effective Lagrangian ( p6D comprising a complete basis of 
scalar dimension-four operators. As a result, Eq. (^) relates Ci to basic, gauge-invariant 
quantities of QCD in a simple way, and is manifestly valid to all orders. 

The formalism developed for QCD in Sections |] and ^naturally carries over to QED. In 
particular, it allows us to evaluate the two-loop 0{aas) and three-loop 0{dal) corrections 
to the decoupling constant (g^ for the fine-structure constant renormalized according 
to the MS scheme. The results formulated in terms of mhin) and may be found in 
Eqs. ( |50| ) and (|5T|), respectively. An appropriately modified LET relates the coefficient 
function Ci^ of the heavy-quark effective ■j'jH coupling to (g-y. Similarly to the case of 
the ggH coupling, this enables us to calculate Ci^ through four loops, i.e. to with 
n = 0, ... ,3. The final results written with rrihlfi) and Mh are listed in Eqs. (^) and 
(p6|), respectively. 



Note added 

After the submission of this manuscript, a preprint ||3^ has appeared in which a^^\Mz) 
is related to a^^^M^-) with four-loop evolution 0] and three-loop matching at the 
charm- and bottom-quark thresholds. In particular, the theoretical error due to the 
evolution procedure is carefully estimated. In Ref. |39|, also the logarithmic terms of the 
three-loop decoupling constant (g, which was originally found in Ref. ||^, are confirmed 
using standard RG techniques. These techniques were also employed in Ref. in order 
to check the logarithmic part of the full result, which was obtained there through explicit 
diagrammatical calculation. 
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Appendix 



A Results for the gauge group SU (Nc) 

In the following, we list the decoupling constants and Q appropriate for the general 
gauge group SU(A^c)- The results read 
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For iVc = 3, we recover Eqs. and (|2|). 



B Decoupling relations for the quark and gluon fields 

In analogy to Eqs. (0) and (Q), the renormalized decoupling constants (2 and Cs for the 
quark and gluon fields, respectively, arise from the relations 



Of course, (2 and (3 are both gauge dependent. Restricting ourselves to the case 
we find in the covariant gauge (H) 
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